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Abstract 

Following the approach of Julien Lesgourgues [1], we analyze the 
mathematical structure of the time co-ordinate of present day 
cosmological models, where these models include a cosmological 
constant term to account for the observed acceleration of the uni- 
verse [2], [3]: we find that in all cases, except for a set of measure 
zero in the parameter space, the time is given by an (abelian) in- 
tegral on a torus; the imaginary period of this integral then gives 
a natural periodicity in imaginary time for the universe; follow- 
ing Stephen Hawking [4], [5], this periodicity may be interpreted 
either as giving a fundamental mass scale for the universe, or (us- 
ing Planck's constant) a fundamental temperature, or both. The 
precise structure that emerges suggests that the structure of time 
can be regarded as an order parameter arising perhaps in a phase 
transition in the early universe; one might hope that this structure 
would be predictable in any fundamental theory. 



Introduction 



The study of cosmology has been transformed recently by the discovery that the 
universe is apparently accelerating [2], [3]. If this is true, the standard cosmo- 
logical models have to be updated to include the acceleration. Typically these 
refined models involve the introduction of a cosmological constant term in the 
formula governing the Hubble function (see below). It is this case that we analyze 
here, although much of the structure we discuss would go through with appropri- 
ate modifications for other kinds of Hubble functions. The space-time metrics in 
question can all be written in the form: 

Here g,^ is a three-dimensional metric of constant scalar curvature k (which is 
naturally compact if «; > and could be either compact or non-compact, for the 
other values of k,). Note that the time co-ordinate t is canonical in that dt is of 
constant size and perpendicular to the surfaces of homogeneity, which are given 
by the equation t — constant. The function a{t) > is the scale factor for the t = 
constant three-geometries. In turn a(t) is related to the Hubble function H{t) by 
the differential formula: da{t) = a{t)H{t)dt. As discussed nicely by Lesgourgues 
[1], the Hubble function H{t) is related to the matter and energy content of the 
universe by the following formula: 

3c^H\t) = SngipR + Pm) - 3Kc'^a-\t) + Ac^. 

Here A is the cosmological constant term, c is the speed of light and Q is New- 
ton's gravitational constant. Also Pr > represents the radiation content of the 
universe and pm > the matter content. For the purpose of this work, we take as 
given that these have the standard characteristic decay rates that p — pRa{tY and 
a — pMO^{t) are positive constants. Then our formula for the the Hubble function 
can be re-written: 

?,(?a\t)H''{t) = 87r^(p + a{t)a) - ^Kc^a'it) + K(?a\t). 

The structure of this relation is clarified if we take a{t) as the basic variable, in 

which case the parameter t is determined implicitly by the differential formula: 

da 

at — -— , where now the Hubble function is regarded as a function of a. In 

aH[a) 

the model given above, H{a) is determined by the formula: 

Sc^a'^H^ = SnQip + aa) - S/^c^a^ + Ac^a^ 
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Using the scale factor a as co-ordinate, the metric now takes the form: 




-9k 



) 



We find it convenient for our exposition to use instead of a, the inverse scale factor 
u — as our co-ordinate. Then the metric becomes: 



Lastly, we absorb a factor of c ^ into H, so we put H — cy, bringing the metric 
to its final form: 



Here both y and u have the dimensions of inverse length. Note that, for an ex- 
panding universe, u decreases in the direction of increasing time. We call y the 
Hubble scale. It is related to the standard time co-ordinate t by the differential 
equation: dt — —{ucy)~^du. Also of interest is the conformal time parameter r, 
defined by the relation dr — —{cy)~^du. It is the time co-ordinate such that dr 
is a conformal Killing vector field, perpendicular to the surfaces of homogeneity. 
Using conformal time, the metric is written g = u~'^{c^dT'^ — g,^), so the metric 
c^dr'^ — g^ represents the conformal structure of the physical metric. We regard 
this conformal structure as fundamental, defining the basic causal structure of the 
universe, the actual metric being given by the conformal factor u'"^ applied to this 
conformal metric. 

We call the right-hand side of the polynomial equation relating y and u the cosmic 
quartic. The equation itself is of a standard type, occurring in the basic theory of 
Riemann surfaces, defining (if u and y are allowed to be complex numbers) a torus 
in the complex (m, ?/)-plane, provided that the quartic has distinct roots, which is 
the generic case. It is the analysis of this torus that will occupy us in this work. 




3kc%^ + Ac^. 
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• In section one below, for completeness, we develop the connection and cur- 
vature of the metric. Our main purpose here is to emphasize that the Hubble 
function is co-ordinate independent and directly measurable, in principle. 

• In section two, we show how to redefine the various variables to formulate 
the cosmic balance equation, essentially identical to that of Lesgourgues 
[1], except that we arrange the definitions such that the key quantities are 
all positive. Whenever we make variable changes, we relate these directly 
to the parameters of the balance equation. 

• In section three we give an overview of the structure of the cosmic quartic. 
We make rescalings of the various variables involved, so that the quartic is 

reduced to the convenient form x'^ + — -— H . Here x — pa~^u. 

32 _ 512 ^ 

Also the parameter is a positive multiple of k pa ^ and the parameter S is 
a positive multiple of Ap^a"^. So we have a two-parameter {R, 5')-space of 
quartics to examine, with the axis R = corresponding to Euclidean topol- 
ogy and the axis S = corresponding to the case of vanishing cosmological 
constant. 

• In section four we show how to parametrize the roots in an efficient way 
and give examples of the various root types that occur. 

• In section five, we bring out the invariant structure of the quartic: it is based 
on using a two-component spinor n^, appropriate ratios of whose compo- 
nents determine the parameter u. We find that the conformal structure of 
the metric takes a natural form: {y^^i^Ad'KBY — G^, where y^^ — —y^^ 
and y^^y'-^^TCATCcH^j^irETTF = 4. Here H^^ is a constant spinor, which 
is symmetric in the pairs AB and CD and is trace-free. The spinor H^j^ 
incorporates the parameters defining the universe. The conformal time is 
then defined by the formula cdr = y'^^TiAdnB. The various special cases, 
where the quartic has repeated roots are treated one by one. 

• In section six, we carefully analyze the general case where the quartic is 
non-degenerate. The parameter r now has two periods, one real and one 
purely imaginary. It is shown how to conveniently make reparametrizations 
to calculate these periods explicitly. 

• In section seven, the results are summarized and the structure of the time pa- 
rameter is discussed. The measurement problem is analyzed and discussed 
in relation to present knowledge. 
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1 General cosmological metrics 



The cosmological structure of the universe is apparently consistent with being 
conformally flat. As described in the introduction, current cosmological models 
use the following metric, which is spatially homogeneous and isotropic: 

du (S> du dr <S> dr 



- r^de ®de -r^ sm\e)d(j) (g) dcj)) 



y'^(u) 1 — Acr^ 

Here u is positive and the Hubble scale y{u) is a smooth positive function, to be 
determined by the phenomenology. In this section we derive the basic properties 

of this metric, for general y{u). We take the domain of the function y to be the 
open interval {Uq, Ui), where < Uq < Ui < oo. Physically u determines the in- 
verse of the radius of curvature of the level surfaces u = constant, so U^^ and 
would then determine the maximum and minimum possible radii, respectively, if 
either or both exist. The co-ordinate r has the range (0, oo), unless k > 0, in 
which case the range of r is the open interval (0, The co-ordinate 9 has the 
range (0, vr) and the co-ordinate cf) has the range (— vr, tt). Also k is an unknown 
constant, determining the spatial geometry: 

• K — gives Euclidean spatial sections. 

• K > gives spatial sections that are three-spheres. 

• K, < gives spatial sections that are hyperbolic three-spaces of constant 
negative curvature. 

Experimentally, it appears that the value k = is presently favored. By an ap- 
propriate constant rescaling the r co-ordinate and the function y, if necessary, it is 
sufficient to restrict to the cases k = 1, k = 0, or k, = — 1. Also, in each case the 
spatial sections could be quotiented by a discrete symmetry group to give a spatial 
section with non-trivial fundamental group. 

The metric g is clearly conformally flat, so has vanishing Weyl tensor. Using 
Maple, the covariant Ricci tensor of g, Rab, takes the following form: 

Rab = i-uyy' + + 2ku^)gab - 2TaTb(yuy' + ku^). 

Here T — —{uy)'^du is a unit timelike future pointing co-vector perpendicular 
to the surfaces of homogeneity u — constant. The Ricci scalar R is: 

R^Q>{-uyy' + 2y'^ + KU^). 
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The Einstein tensor Gab is then: 

R 

Gab = Rab - i^Qab = {'^uyy' - - Ku^)gab - 2TaTb{yuy' + ku^). 

• The metric is flat if and only if Gab vanishes, which gives the relation: 

2 2 

y — —Ku . 

This entails that k < 0. Then y = u\f^^ and is defined for all positive u. 

• The metric is of constant curvature, if and only if Gab is proportional to gab, 
which gives the relation: 

y''^C-KU^. 

Here C is an integration constant. When y^ — C — kv?, the Einstein tensor 
reduces to Gab = -^Cgab- 

- If K = 0, we need C > and then y is defined for all positive real u. 

- If K > 0, we need C > and then y is defined for the li-interval 

- If K < and C > 0, then y is defined for the li-interval (0, oo). 

- If K < and C < 0, then y is defined for the li-interval (VCkt^, oo). 

Henceforth we assume that the metric is not of constant curvature on any open set, 
so that the quantity yuy' + ku^ is never zero on any open set. 

• Taking the divergence of both sides of the equation for the Einstein tensor, 
we get, using the fact that daU — {du)a — —uyTa and the fact that the 
Einstein tensor has zero divergence: 

= -2{yuy' + KU^){nd^Ta + Tad'^n) + ?>uyn{KU^ + y^)' 

= -2{yuy' + Ku'){WTa + Tad^T, - SyT,), 
= T,9"T„ + T„9"T,-3|/T,. 
Taking the trace of this equation with T'^, we get the relations: 

d'^Ta = 3y, 

Ta&'n = 0. 

This tells us that is a parallelly propagated unit geodesic timelike vector field. 
Also measurement of the quantity |9aT", one third of the co variant divergence of 
the vector field, T", determines the Hubble scale factor y. 
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2 The cosmic quartic and the balance equation 



The specific metric discussed in cosmology takes the fomi: 

3cW dr^ 



ds — u 



r^{d9^ + sin" {9)d<j)^) 



_8ng{pu^ + (7^3) - 3kc%2 + Ac2 1 _ 
So we have the following fonnula for the Hubble scale y: 

3c^y^ = SnQipu^ + au^) - 3nc\'^ + Acl 

Here c > is the speed of light, ^ > is Newton's gravitational constant, p > 
measures the radiation density and a > measures the matter density. Finally A is 
the cosmological constant term (which, if positive, dominates the geometry when 
u is small). For the remainder of this section we will assume the generic case that 
A 7^ 0. We can rewrite the equation for y in terms of present day quantities: 

3cY = SnQpo f — 1 + SnGao ( - 3kc%1 ( + Ac^ 

\UoJ \UoJ \UoJ 

= F{X), F{X) = flRX"^ + QmX^ - Kfl^X^ + eflA- 

Here — k, — 1, uq is the present value of u and yo is the present day value 
of y. Also we have: 

Z=^, X=-, y^ = J^. «2^Mk, 
P = PoUo ^ o n 3 = o IAI09 ' (^ = (^oUo - - 



This last equation is called the cosmic balance equation, reflecting the relative 
contributions of matter (0,m), radiation (Qr), topology and the cosmological 
constant (JIa)- The quantities Z, X, Qr, Qm, and fl^ are all dimensionless and 
positive. 
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3 A unified treatment of the cosmic quartic 

Consider the equation for the Hubble parameter: 

Make a change of variables: 

xa 

• u — — , 

P 

^ 64AcV^ 

Note that R and S are positive multiples of k, and A, respectively and x, Y, R and 
5" are pure numbers (dimensionless). When k, and A are non-zero, we have: 



u — 



x\A\ 2 / Hm^k 



F|A|5 



3 1 



R 



32/s: / n^R^K 



9 V 
512A /Q|Qa 

In terms of these variables the Hubble parameter equation becomes: 

Q/? ^ 

3F2 = Eix) = x" + - — a;^ + — . 

^ ^ > 32 512 

From the present viewpoint, to precisely determine the metric of space-time, the 
very difficult experimental task is to measure accurately the pure numbers R and 
S. 
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Note that a priori the parameters R and S are continuous variables. In particu- 
lar, if the sign of R is determined, then we will know the topological type of the 
spatial geometry. Currently it seems to be believed that the universe is spatially 
Euclidean. However it is rather dangerous to take this seriously, as this requires 
that R be exactly zero. A small deviation in R from zero, of either sign would 
drastically alter the spatial character of the metric, without substantially altering 
its temporal evolution on time- scales comparable with the present age of the uni- 
verse. In this section, we give a unified treatment of the cosmic quartic H{x), with 
no prior assumptions about the values of the quantities R and S. First consider the 
behaviour of the roots of H as the parameters R and S are varied. The resultant 
J{R, S) of H and its derivative with respect to x factors as: 

J{R, S) = 2-^^S{6AS^ - A32S{3R^ + 12R + 8) + 729R%9R + 8R^)). 

Solving the equation J{R, S) — 0, for S in terms of R, gives three roots, {0, S+, S-}, 
where S± are given as follows: 



Deleting the zero set of J divides its complement in the {R, S) plane into five con- 
nected open regions, where the character of the roots is the same in each connected 
component. We discuss each in turn: 

• Region I: < S < S.,R> 0. 

In this region H has four real distinct roots, two positive and two negative. 

• Region H: S_ < S < S+ and S > 0. 

In this region H has two real distinct roots, both negative and two non-real 
roots. 

• Region III: S > S+ and S > 0. 

In this region H has all its roots non-real. 

• Region IV: < with R < -lor R> 0;or S+ < S <0, with -1 < R < 
-|; or 5 < S_, with -1 < i? < 0. 

In this region H has one real root of each sign and two non-real roots. 

• Region Y: S- < S < and -1< R < 0. 

In this region H has one positive real root and three distinct negative roots. 
This is the only bounded region. 




9 



The edges of the regions correspond to cases that the roots are repeated. 

• The origin R — S — is a boundary point of all regions except for region 
ni. The roots are 0, 0, and —1. This is the Euclidean case with A = 0. 

• All the regions except for region I have a common boundary point with 
R——^ and S — 0. The roots in this case are — | and 0, each double roots. 

• Regions IV and V have a boundary cusp point at i? = — 1 and S — 
The roots in this case are a simple root of | and a triple root of — |. 

Deleting these three special points, the boundaries split into seven open connected 
curve components: 

• I-II: The curve separating regions I and II is given by the equation S = S-, 
with R > 0. Here the roots are all real with two equal positive roots and 
two unequal negative roots. 

• II-III: The curve separating regions II and III is given by the equation S = 
S+, with i? > — |. Here the repeated roots are negative and there are two 
non-real roots. 

• III-IV: The curve separating regions III and IV is given by the equation 
S = 0, with i? < — |. Here the repeated roots are zero and there are two 
non-real roots. 

• IV- VA: The first curve separating regions IV and V is given by the equation 
S — S-i-, with — 1 < i? < — |. Here the repeated roots are negative and 
there are two other real roots one of each sign. The repeated roots are the 
most negative. 

• IV- VB: The second curve separating regions IV and V is given by the equa- 
tion S — S^, with — 1 < i? < 0. Here the repeated roots are negative and 
there are two other real roots one of each sign. The non-repeated negative 
root is the most negative. 

• II- V: The curve separating regions II and V is given by the equation S = 0, 
with — I < i? < 0. Here the repeated roots are zero and there are two other 
real roots, both negative. 

• I- IV: The curve separating regions I and IV is given by the equation S* = 0, 
with R > 0. Here the repeated roots are zero and there are two other real 
roots, one of each sign. 
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We now summarize the structure of the quartic according to the signs of A and n. 
We first look at the cases where one or both of these quantities vanishes. 



In the case that A = 0, we have S — Q. Then we have: 



H{x) ^x^{^'^ + x- , with roots |o, 0, ^(-4 ± Vl8i? + 16)| 

- The Euclidean point is then R — Q and H{x) — x^{x 

- The elliptic case is > 0, so is the I-IV boundary. Then H{x) has a 
positive and a negative real root. 

- The hyperboUc case has i? < 0, so has three parts: 

* the 11- V boundary, with — | < i? < 0, so has two negative 
real roots. 

* the special point R = — |, when H{x) = \x'^{2x + 1)^. 

* The 111-lV boundary, with R < so H{x) has two non-real 
roots. 

• The Euclidean case with A > has i? = and breaks up into three sub- 
cases: 

- If < (S" < 54, then we are in the interior of region II. 

- \f S > 54, then we are in the interior of region IE. 

- If S — 54, then we lie on the II-III boundary curve. In this case we 



have: 



27 1 

H{x) = + + — = — (16a;^ - 8a; + 3)(4a; + 3)^ 

ZOO ^00 

^ {Ax + 2>f{Ax-l+iV2){Ax-l-i\/2). 



256 

• The Euclidean case with A < has i? = 0. The quartic then lies in region 
rv, so has one real root of each sign and two non-real roots. 
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Finally we summarize the various generic cases where neither A nor k vanishes. 

• The elliptic case with A > corresponds to the open first quadrant of the 
{R, 5')-plane. Here we have all of region I and the parts of the I-II boundary 
and of region n in the first quadrant. 

• The elliptic case with A < corresponds to the open fourth quadrant of the 
{R, yS) -plane. Here we have the parts of region IV in the fourth quadrant. 

• The hyperbolic case with A > corresponds to the open second quadrant 
of the (i?, 5) -plane. Here we have the parts of regions II and HI in the sec- 
ond quadrant and the parts of the II-III boundary curve in the open second 
quadrant. 

• The hyperbolic case with A < corresponds to the open third quadrant of 
the (i?, 5')-plane. Here we have all of region V, both boundary curves IV- 
VA and IV-VB, the cusp point = and the part of region IV 
open fourth quadrant. 
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4 Parametrizing the roots of the cosmic quartic 



For many purposes it is desirable to have good control over the roots of the quartic. 
One method to do this is as follows. Given R and S, introduce a new parameter 
P 7^ by the implicit formula: 

S^^{P-l){9R-SP + 8f. 

Multiplying out, this gives a cubic equation for P, which is straightforward to 
solve for P in terms of R and S by standard methods: 

= g{P) = 64P3 - 48P\3R + 4) + P{3{3R + 8){9R + 8) - 8S) - {9R + 8)^. 

Note that for given R and S, we have g{0) = -{9R + 8f and g{l) = -85", so, if 
if (S" > 0, the equation g{P) = has a solution P > 1, whereas if jS" < 0, there is 
a solution in the interval < P < 1, unless P = — |. 

• In the exceptional case that P = — |, we have a direct factorization of the 
quartic: 

H(x) = x'^+x^^+\x^+^ = -^(32x'^+16x+iV2S)(32x'^+16x-iV2S). 
^ ^ 4 512 1024^ ^ 

The roots are {^(-2 ± \[7+2d\/^), ^(-2 ± \l A - 2i\/2S)} . 

- If 5" > all roots are non-real, with equal sizes for their imaginary 
parts and with two positive real parts two negative. 

- If yS" = 0, the roots are and — \, each a double root. 

- If —2 < 5" < 0, then all four roots are real, with three negative and 
one positive. 

- If 5" = —2, then there is a double root of —\, a third more negative 
root — |(v^ + 1) and a positive root |(V^ — 1). 

- If S" < — 2, there are two real roots, one of each sign, and two non-real 
roots each with real part — |. 

Henceforth we assume the generic case that P 7^ — |. The cubic for P simplifies 
if we write S = ^(T — P^), where T is real. Then the critical points of g{P) are 

|(6P + 8±3Vr). 
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• If < -R^, then T < 0, so g{P) has no real critical points. Then g{P) is 
everywhere increasing and its real root is unique, is larger than 1, and varies 
smoothly with R and S. 

• If 5" = —"^R^, then T = and g{P) has one critical point only, which is 
also its inflection point, so g{P) is everywhere increasing, so its real root is 
unique. Explicitly, we have, for this case: 

P = 1 + ^(1- + 

So for R > 0, P lies in the range | < P < 1 and decreases towards the 
value j,asR gets large. 

• If 5" > -f i?^ then T > R^. Then the value of g{P) at the largest critical 
point P = |(6i? + 8 + 3VT) is f (i? - VT)(T + A{R + 1){R + ^/T)). 
Since T > R^, this value is plainly negative. So in this case g{P) has a 
unique real root, which is larger than ^(18i? + 24 + VSl-R^ + 245*) and 
which varies smoothly with R and S. 

The resultant of g{P) with its derivative factors as —2^^J{R, S), so the roots of 
g{P) change character when those of H{x) do. 

• At (0, 0) in the {R, 5')-plane, all three roots are 1. 

• At (— |, 0) in the (i?, 5')-plane, one root is 1 and the other root is repeated 
and is P 0. 

• At (-1, -^) in the (i?, 5) -plane, aU three roots are P = |. 

• In region I, there are three real roots for P, all positive and larger than 1. 

• In region II, there is one real root for P, which is positive and larger than 1. 

• In region 111, there are three real roots for P, with exactly one positive; this 
root is larger than 1; except along the vertical line P = — |, when one root 
is zero, one is positive and larger than one and one is negative. 

• In region IV, there is exactly one real root for P, which is positive and 
smaller than 1, except along the vertical line P = — |, when one root is 
zero and the other two are non-real with real part |. 
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• In region V, there are three positive real roots for P, all smaller than 1, 
except along the vertical line R — when one root is zero and the other 
two are real and smaller than 1. 

• On the I-II boundary curve there are two equal real roots and a third larger 
real root, all positive. 

• On the Il-ni boundary curve there are two equal negative real roots and a 
third positive real root, which is larger than 1. 

• On the ni-IV boundary curve there are two equal negative real roots P — 
l(9R + 8), with < -| and the third root is P = 1. 

• On the boundary curve IV- VA, all three roots are positive and smaller than 
1, with the two smallest equal. 

• On the boundary curve IV- VB, all three roots are positive and smaller than 
1, with the two largest equal; except at the point (— |, —2) of this curve, 
when one root is P = and the other two are each P = |. 

• On the II-V boundary curve there are two equal positive real roots P = 
|(9P + 8), with -| < P < 0, both smaller than 1 and the third root is 
P = 1. 

• On the I-IV boundary curve there are two equal positive real roots P = 
|(9P + 8), with P > 0, both larger than 1 and the third root is P = 1. 

Replacing S by its expression in terms of P and R, the cosmic quartic factorizes 
readily: 

4 , 3 ^-^^2 , S _ Q+Q- 



Hix) = x^ + x-^ - — 



32 512 4096 ' 
g± = Ux^ + 32x - 9P - 8 + 8P ± (Vp)-^(8(4x + 1)P - 9P - 8). 
The roots are then: 

i (^-2 - 2Vp±^ (1 + (v^)-i)(3^pTT+ 1 - 2Vp)(3v/pTT - 1 + 2VP)^ , 
^ (^-2 + 2Vp± ^J{1- (/P)-i)(3v/pTT + 1 + 2VP){3VrTi - 1 - 2VP)^ . 
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We give at least one example for each case: 

• At (0, 0), H{x) = + = x-^{x + 1), with roots {0, -1}. 

• At(-|,0),ii'(x) = |x2(2x+ 1)2, with roots {0,-|}. 

• At(-1, f),//(x) = 2-i4(128xVl92x-9)(128x2-64x-3),withroots 

l8' 8J' 

• At the point (3, §) in region I, H{x) = 2-^^(128x2 + 192x - 9)(128x2 - 
64x - 3), with roots {^(4 ± a/22), ^(-4 ± 3V2)}. 

• At the point (3, ^) in region 11, 

H{x) = 2-^3-2(96x2 - 96x + 37) (48x2 + 96x + 37), 
with roots {^(-12 ± V^), ^(12 ± iV78)}. }. 

• At the point (0, ^) in region 11, 

H{x) = 2-103-2(96x2 + 120x + 25) (96x2 _ 24^^ + 25), 
with roots {^(-3 ± V3), ^(3 ± iV^l)}. 

• At the point (-|, ffH^) in region 11, 

H{x) = 2-1^3-^(4608x2 + 4896x + 1207) (4608x2 - 288x + 71), 
with roots {^(3 ± iy/US), ^(-51 ± i^l87)}. 

• At the point (3, 3267) in region 111, 

H{x) = 2-9(16x2 - 32x + 33)(32x2 + 96x + 99), 
with roots {|(-6 + ±3iV6), |(4 + ±iVri)}. }. 

• At the point (0, ^) in region 111, 

H{x) = 2-10(32x2 - 48x + 45) (32x2 + 80x + 75), 
with roots {§(-2 + ±iV2), |(2 + ±iV6)}. 

• At the point (-|, ^|||^) in region 111, 

H{x) = 2-1^(512x2 - 256x + 123)(512x2 + 768x + 369), 
with roots {§(-8 ± 3iV2), ^(8 ± iVTS2)}. 

• At the point (3, -^) in region IV, 

H{x) = 2-12(64x2 + 16x + 33) (64x2 + 48x - 99), 
with roots {§(-1 ± 2v^), |(-1 ±4iv^)}. 
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At the point (0, - y) in region IV, 

H{x) = 2-i°(32j;2 + 8x + 3)(32x^ + 24a; - 9), 
with roots {§(-1 ± V3), ± iVb)}. 

At the point (— 1, — ^) in region IV, 

H{x) = 2-14(128x2 + 48a; - 9)(128a;2 + 80a; + 15), 

with roots {^(-1 ± \/3), ^(-5 ± iy/b)}. 

At the point (— |, — |) in region V, 

H{x) = 2-i°(32a;2 + 24x + 3)(32a;2 + 8a; - 1), 
with roots {i(-l ± \/3), i(-3 ± \/3)}. 

At the point (3, ^) on the I-II boundary curve, 

H{x) = 2-i2(64a;2 + 112a; + 21)(8a; - 3)^, with roots {§, |(-7 ± 2^7)}. 
At the point (3, ^) on the 11-111 boundary curve, 

H(x) = 2-12(64^2 _ + 45)(8^ + 9)2^ with roots {-|, |(5 ± 2iV5)}. 

At the point (0, 54) on the II-III boundary curve, 

H{x) = 2-«(16a;2 - 8a; + 3) (4a; + 3)^, with roots {-|, |(1 ± iV2)}. 

At the point (— |, fH) on the II-III boundary curve, 

H{x) = 2-i6(256a;2 - 32a; + 9) (16,t + 9)^, with roots {-^, ^(1 ± 2i^/2)}. 

At the point (— y|, ~|^) on the IV- VA boundary curve, 

H{x) = 2-20(32a; + 5)(32a; - 3)(32a; + 15)^, with roots {-|§, |)}. 

At the point (-^f, -§f|) on the IV- VB boundary curve, 

H{x) = 2-20(1024a;2 + 448a; - 63) (32a; + 9)^, 
with roots {-^,^(-7±4\/7)}. 

At the point (3, 0) on the I- IV boundary curve, 

H{x) = 2-V(32a;2 + 32a; - 27), with roots {0, |(-4 ± VTO)}. 

At the point (— ^, 0) on the II-IV boundary curve, 

H{x) = 2-V(512a;2 + 512a; + 63), with roots {0, ^(-16 ± v^)}. 

At the point (—1, 0) on the lU-IV boundary curve, 

H{x) = 2-5a;2(32a;2 + 32a; + 9), with roots {0, |(-4 ± iV^)}. 
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5 The invariant approach to the cosmic quartic 



We introduce a (real) two-component spinor tta, the ratio of whose components 
represent the variable x. Let H^^'^^ be a fixed spinor, obeying the following 
symmetry conditions: 



TT CD _ TT CD _ TT (CD) 

^AB — -"(AB) — ^AB ' 



H 



AB 



AC 



0. 



Introduce a skew spinor symplectic form y^s = —Uba and its (skew) inverse y"^^, 
defined in terms of Hj^^^^ and tte by the relations: 

TT EF AB o 

VABVCD ^ T^[CT^[AnB\D\ '^E'^F, UABV = ^, 

y^^y^^n^ncH^^^^nEn^ = 4. 

Note that yAB{knc) = k^yAB{T^c) and y'^^{k'Kc) = k'^y'^^iTic). Then the 
conformal metric G of cosmology achieves the following canonical form: 



G = (y^^nAdTTBy 



r2(d^' + sin'(^)#'), 



Here all the variables tt^, r, 6, (p, y^^ and Hj^^'^^ are dimensionless. Also — 
K. The information of the cosmic quartic is encoded in the constant spinor i^^^*^^. 
In turn the spinor H^^'-'^ may be regarded as tensor Hj", where the Latin indices 
a <-> (AB) are three-dimensional, such that = and Hab — H(^ab), so Hj" 
has five degrees of freedom. The metric G is invariant under the natural action of 
GL(2, M) on all the spinors. Since the spinor iJ^^'"^ has equal numbers of upper 
and lower indices, the action of GL(2, M) factors through the action of SL(2, M) 
and even of PL(2, M), or the pseudo-orthogonal group, SO(l, 2). The group action 
has three degrees of freedom, giving the orbit structure two degrees of freedom. 
With respect to a suitable basis, the matrices Gq for the covariant SO(l, 2) metric 
and H for the cosmic quartic take the following forms: 



Go^ 



1 
0-2 

1 



H 




_s_ 

512 



-1/2 

32""' 



4-. 



1 

1 

4 



It is easily checked that H is trace-free and GqH is symmetric, as required. 
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In this language, a spinor with components (p, q) is encoded in a null vector 
s° with components {p'^,pq, q^). Also a co-spinor tt^ with components y) is 
encoded in a null co-vector pa with components (x^, 2xy, y"^). Note that we have 
s"'Pa = {cf'^t^aY- The characteristic polynomial p( A) of the matrix H is: 

The condition that the roots be distinct is precisely the vanishing of the polynomial 
J{R, S). Also in this language, we have ha proportional to the co-spinors (x, 1) 
and {up, a). So if we take as basis spinors, the spinors = (1, 0) and {3^ — 
(0, 1), then we have: 

ttbch^ pu 

The (dimensionless) conformal time parameter is defined by the formula: 
The (dimensionless) time parameter is defined by the formula: 

A 

-1 TT^Q; ^r) 

dt ^ X dr ^ — -^y"^ TTcdiTD- 

The philosophy now is that the equation y"^^ y^ t[ at^cH ^d^^ et^ f = 4 generi- 
cally defines a torus (i.e. provided the tensor has distinct eigen-values), when 
y"^^ and tta are allowed to be complex. Then there are two independent integrals 
defining each of the quantities t and r (i.e. four integrals in all), which can be 
chosen so that, for each quantity, one of the integrals is real and the other com- 
plex. The periodicity in the variable r is the natural geometric periodicity and 
characterizes the complexified conformal geometry; for example all the Green's 
functions of solutions of conformally invariant hyperbolic wave equations, such 
as Maxwell's equations, on the cosmological background will exhibit this period- 
icity. Then the corresponding integrals for the co-ordinate t will give the lifetime 
of the universe for the real integral (if it is finite) and a natural periodicity in imag- 
inary time for the imaginary integral. Note that the lifetime is always finite in the 
past going to the big bang as — > oo. 
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Temporarily going back to our initial variables, we have: 



t = 



^/Scd' 



u 



y/Scdu 



^/Sngifm'^ + au^) - Skc^u'^ + Ac^ 

Note that if r is used as a parameter and if then the variable u is known as a 
function of r, the integral for t simplifies to: 



t 



_ f dr 

~ J wy 



Passing to the dimensionless versions of these relations, the quartic is now written 
in terms of the homogeneous variables x and y as: 

H{x, y)^x' + x'y - ^x'y'R + ^v'S = //"Wo- 

Then the dimensionless times r and t are given by the following scale invariant 
integrals: 

y{xdy - ydx) 



t 



T — 



X\l x'^ + x^y — ■^x'^y^R + 5^2^"''^ 



{xdy — ydx) 



^x^ + x^y-^x^ym+^y^S 

In the generic case, we find that u varies periodically as a function of the real 
variable r oscillating between two simple real roots of the quartic. Again if we 
parametrize the ratio x : y in terms of r, we may rewrite the integral for t as: 



t 



ydr 



X 



Note that from the point of view of complex integrals, the integral for t always 
has an additional pole at x = 0. This introduces an ambiguity in the real integral 
if S" < and in the imaginary integral if S" > 0, corresponding to 27r divided by 
the modulus of the square root of the quartic at (x, y) — (0, 1). 
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We analyze the behavior of the conforaial geometry according to the root structure 
of the cosmic quartic, where we write the spinor components of vr^ as (x, y) and 
the quartic as H{x,y). For completeness, we include a few cases that wouldnot 
be realized physically. 

• In the (unphysical) case that the quartic is a perfect fourth power, H{x, y) — 
{ax + hyY, with a 7^ 0, we have the parametrization: 

[x,y\ = [1 - afer,aV]. 

The conformal time r ranges over the whole real line. Then we have: 

t^h-'^{-ahr -\vL{\ahr 

If a = 0, so H{x, y) = c^y^, with c > 0, we have the parametrization: 
[x, y] = [rc, 1] and then t — —c"^ ln(|r|). Again the conformal time ranges 
over the whole real line. In each case, as r ^ ±00, the root is approached. 
Also, in each case, all ratios x : y with the exception of the triple root are 
traced out once, as r varies. 

• In the case that the quartic has a triple root (the two physical possibilities 
are the origin, {R, S) = (0, 0) and the cusp point, {R, S) = (-1, -y)), 
there are two types of quartics: 

- First we may take H{x, y) — x^{x + 2ay), where a 7^ (the physical 
case is a = |). Then we may parametrize (x, y) as: 

[x,y] = [2a,aV2-l]. 

Here, a priori, the variable r ranges over the whole real line. Note that 
X + 2ay = la^T^. As r ^ ±00, the triple root is approached. At 
T = the simple root is reached, but not crossed. 
We need ^ to be positive, so if a > 0, we need a^r^ > 1 and if a < 0, 
we need a^r^ < 1. Then we have: 

The physical case has a > and | > 0, so occupies the interval r > ^ 
and the age of the universe is infinite. 



21 



- Second, we may write the quartic as: H{x, y) = {x — 3ay)^{x + ay), 
where a 7^ (the physical case is a = — |). Then we may parametrize 

{x, y) as: 

[x,y] = [a(12aV + l),4aV - 1]. 

Here, a priori, r ranges over the whole real line. We have x — Say = 4a 
and X + ay — 16a^T^. As r — > ±00, the triple root is approached. At 
r = the simple root is reached, but not crossed. 
We then have: 

t = ^^(3aT — 2-\/3arctan(2a-\/3T)). 

For the physical case that a < 0, r may be taken to range between 
^ and — so the physical age T of the universe is finite and is (in 



dimensionless units): 



T=^(4.v^-9). 



In particular for the case a = — |, so {R, S) — (—1, — y), we have 

2567r\/3 

Thus, in each case, the ratio x : y for the conformal geometry moves from 
the triple root as r ^ —00 to the single root at r = and then retraces its 
steps, going back to the triple root as r 00. However the evolution in 
physical time t occupies only a part of this full trajectory. 

In the case that the quartic is a perfect square, with real roots, (the physical 
case is {R, S) = (— |, 0)), we may write the quartic as follows: 

H{x, y) — x^{x + ayY, where a 7^ 0. 

The physical case is a = |. Then we may parametrize {x, y) as: 

(x,y) = [a,e--l]. 

We have x-\-ay — ae"^. Then, as r — > —00, one of the roots is approached 
and as r — > 00, the other is approached. Note that in this case, there is a 
natural periodicity in imaginary r, with period 2niar^. In particular for the 
physical case that H{x, y) = x^{x + |)^, we have period Airi. Physically, if 
we take a > 0, we need m > 0, so r > 0. Then we have: 

t^a-^e"^ -ar). 

So t has infinite range. 
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• Alternatively, again for the case that the quartic is a perfect square, we may 
write the quartic in a way that allows the roots to be complex: 

H{x, y) = (x^ - ec^y^Y, where c> 0, = 1. 

- If e = 1, then we may parametrize as follows: 

[x,y] — [ccosh(cT), sinh(cT)]. 

We have x + cy = ce^'^ and x — cy = ce^'^'^. Then, as r — — oo, one of 
the roots is approached and as r ^ oo, the other is approached. Note 
that in this case, there is a natural periodicity in imaginary r, with 
period Tric~^. This case is equivalent to the previous case, modulo 
an appropriate SL(2, M) transformation. Under that transformation, 
we find that c = ^, so the periods are in agreement. Also we have 
t — ln(cosh(cr)). Here r > 0, so the universe has infinite age. 

- If e = —1, then we may parametrize as follows: 

[x,y] — [ccos(cr), sin(cT)]. 

Then the conformal geometry is periodic, with period ttc"^, covering 
all ratios x : y in the period. This case is not physical. 

• In the case that the quartic has a single repeated root, generically it takes the 
form: 

H{x,y) = ^{x + 2ay)^{bx - {b - l)ay){x + (6- l)ay). 

Here a 7^ and 6 7^ 0. Also, without loss of generality, we may take b to 
be real, with < |6| < 1 and 67^ |, or 6 to be complex with \b\ — 1 and 
9(6) > 0. Putc= |a|y|6-i(3- 6)(36-l)| > 0. 

- We first assume that b is real and < 6 < |. 
We parametrize (x, y) as follows: 

[x, y] = [-2a(6 - 1)(2(1 - b) + z{b + 1)), 1 - 66 + 6^ + z{b'' - 1)]. 
Then the equation for dr reduces to: 

c\dTf = - c= |a|Vl&-H3-&)(l-36)| >0. 

1 — Z'^ 
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The left-hand side is positive, so necessarily \z\ < 1. Then we may 
take z — esm.{cT), where — 1, giving: 

[x,y] = [-2a(6-l)(2(l-6)+esin(cr)(6+l)),l-66+6^+esin(cr)(62-l)]. 

When sin(cr) = e we are at the simple root 2(6 — 3)[a{b — 1), b]. 
When sin(cT) = — e we are at the simple root 2(1 — 36)[a(l — 6), 1]. 
So the conformal metric oscillates forever between the two roots, never 
crossing either, and avoiding the double root completely (since x + 
2ay = 2a(l — 36)(3 — b) ^ 0). Its period is 2nc~^. For the physi- 
cal metrics we may take a — jp^i^- Then we have z in the interval 
(— 1, 6o) where we have: 

Putting z = j^, the time function is then: 

arctan(p)(6 1)^ /, .n4, 3,, ^, , / (b — 3)p 
t = ^^^^ — h 6+1 6-1 -1 arctan , ^ 

V&^(&-3)(36-l) y^^3b-l)ib-3)) 

The time T elapsed during one cycle of the universe, going from the 
big bang to its maximum radius and back again to zero radius is: 



26-i(6+l)^ I (l-6)-iarctan(6-^)+ 



arctan ' ' ^ 



V(&-3)(36-l) 

This diverges as 6 — > 0+, when the double root becomes a triple root. 
- We next assume that 6 is real and |<6<1, or— 1<6<0. 

We parametrize (x, y) as follows: 

[x, y] = [-2a(6-l)(4(l-6)z+(z2+i)(6+i)), 2{l-6b+b'')z+{z^+l){b''-l)] 
Then the equation for dr reduces to: 

c\dTf = ^, c= |a|Vl&-'(3-6)(l-36)| > 0. 

z 



24 



We may take z — ee"", where = 1. So we have: 

{x,y\ = [-2a(6-l)(4(l-6)ee^"+(e^^"+l)(6+l)),2(l-66+6^)ee""+(e2^"+l)(62-l 
Also we have: 

x + 2ay^ 4:az{3b - l)(b - 3), 

{bx - (6 - l)ay) = -a{b - 1)(36 - + l){z - if, 

x + {b- l)ay = a{b - l){b - 3){b + l){z + if. 

As T — oo, we approach the double root. At r = we are at the 
root - = l^-llfi: if e = 1, or we are at the root - = (1 — b)a if e = —1. 

y b y ^ ' 

This root is not crossed and as r — > cxd, we go back to the double root. 
The motion is periodic in imaginary r with period 27rc~^i. Note that 
the values of the ratio ^ at r and — r are equal. With a = j^^^ and 
— 1 < 6 < 0, the physical time function is: 

(b+lf (. „ „ ^ ^ /(2(l-6) + ^(l + 6)) 



2bc{b-l)y ^ ^' 1^ V(3&-l)(&-3) y 

With a = (^^ri)2 and | < 6 < 1, the physical time function is: 

t ^ #±li; U-b) H\A) + 2arctan { W^^±I^\ ^Z^Zb^ 

2bc{b-i) y ^ ^' \^ y 

In either case, the universe has infinite age, since t diverges as 2; — > 0, 
so as the double root is approached. 

- Finally we assume that b is complex with |6| = 1, with Q{b) > 0. 

1 + tp 

Write b = , where p is real and positive. Explicitly, we have 

1 — zp 

i(l - b) 

p — ^ ^ ^ . Then make the parametrization: 

[x,y] = [-2ap{z'^ - Azp - l),pz^ + 2z{2p^ + 1) - p]. 
With this parametrization, the equation for dr becomes: 

c^dT^ = ^, c=\a\yJ\b-\'i-b){l-Zb)\ > 0. 
z 
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So 2; = ee"", where = 1 and r ranges over all the reals. Then we 
have X + 2ay = Aa,z{Ap'^ + !)• So the metric cycles once around, 
starting and ending at the double root as r ±00. As it does so, all 
values of the ratio a; : 1/ are traced once, except fox x : y = — (2a)~^, 
the double root. Note that this metric has a natural period of 27rc~^i 
in imaginary r. In this case, the points parametrized by z and z' agree 
\f zz' = —1, so we can also make a "Mobius" identification in the 
complex, identifying r with —T + 7ric~^. The time function is: 

t = -4p"^(l+/)"^(i+4p2)-| ^pin(|;2|) - 2(1 + 4p^)5arctanh 

Again the physical age of the universe is infinite. 
For these cases, the corresponding values of the parameters R and S are: 
_ [326(3 - 26 + 362) 20486^(6 - l)^" 

^ ' ^ ~ ~ [ 9(6 + 1)4 ' (^^r^p . ■ 

- When 6 ^ I we go to the cusp point (—1, — y). 

- When 6 ^ 1, we go to the point (— |, 0). 

- When 6 — > 0, we go to the point (0, 0). 

- The 6-interval (|, 1) traces out the boundary curve IV- VA in the third 
quadrant. The period goes steadily decreases as R increases, going 
from 00 at the cusp point, to Ani at the point (— |, 0). 

- The 6-interval (0, |) traces out the curve IV- VB in the third quadrant. 

- The 6-interval (—1, 0) traces out the curve I-II in the first quadrant. 
The period decreases from 00 as — > 0+ to zero as i? — > 00. 

- When b is complex, not real and |6| = 1, we trace out the boundary 
curve Il-ni boundary curve. Put 9o — arctan(v^) and denote the 
argument of 6 by ^, so < ^ < tt. The ^-interval (0,^0)5 or < 
P < traces out the then we are in the part of the II-III boundary 
curve in the second quadrant. The ^^-interval {6q,u:), ox p > ^ traces 
out the part of the II-III boundary curve in the first quadrant. When 
9 — 9q, ox p = we are at the Il-III-boundary point on the -S-axis at 
(0, 54). The period steadily decreases as R increases, starting at 47ri 
at (^|, 0), going down to |7ri"\/2 at (0, 54) and going towards zero as 
R ^ 00. 
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In summary, along the upper curve S — S+{R), except where it meets 

the lower curve S = S_{R) at the cusp point, and along the part of the 

lower curve S = S-{R) in the first quadrant, the metrics have a periodic- 

• • • • T u u • J • 2m\b+l\^ 

ity in imagmary r. In each case the penod is — , Note 

^\b{3b-l){b-3)\ 
that this goes to infinity as 6 — > |, which gives the cusp point (i?, -S") = 
(—1, — Y),and as 6 — > 0, which gives the origin. Also as b ^ 1, the period 
goes to Ani. This is agrees with the imaginary period found previously, for 
this case, the case {R, S) = (— |, 0). 

The non-generic cases with two repeated roots are the cases with S — 0. 
For these, we may write the quartic H{x, y) as follows: 

H{x^ y) — x^{x^ — 2bxy — ecy^). 

Here b and c are real, with c > and b^ + ec 0. Note that it! and e have 
the same sign. 

- When e = — 1, (corresponding to i? < 0, which gives the II- V bound- 
ary curve for c < 6^ and the ni-IV boundary curve for c > b^), we 
make the parametrization: 

[x, y] = [Acz, 4(6^ -c) + z^ + Abz]. 

Note that x'^ — 2bxy + cy'^ = c{z'^ — 46^ + 4c)^. Then we have cdr'^ = 
z^'^dz^, so we may take z = ae^^, where = 1. When r — > ±oo, 
the ratio x : y goes to the double root. If 6~c < 0, it simply goes 
through all values of the ratio x : y, except for 0, as r ranges through 
the real numbers. If, instead, 6^ — c > 0, then x : y goes to the root: 
2c : aV^^ — c+ 1, when e^^"^' = 4(5^ — c). So the solution goes from 
the double root, to a single root, which it does not cross, and then back 
to the single root. In both cases, the solution is periodic in imaginary 
T with period 2mc~^. In terms of the physical variable R, the period 
is 3^^2^ • In particular the period goes to infinity as i? — > 0~ and the 
period goes to Airi as —> — |, both as expected. As — > — oo, the 
period goes to zero. The physical time function is: 

^{-z^ - Abz\n{\z\) + Ab"^ - Ac). 

Azc-2 

The physical age of the universe is infinite. 
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- When e = 1, (corresponding to it! > 0, so the I-IV boundary curve), we 
make the parametrization: 



[x, y] = [-C, b + zVc + b'^]. 

Note that x"^ — 2bxy — cy"^ — c{l — z'^){c + b"^). Then we have 
cdr^ = -, so we may take z — cos(tv^). The solution is pe- 

riodic in r with period tq — ^ — j^^- For n any integer, at 
T = 2nro, the solution is at a single root, then goes to the double root 
when cos(ry^) = — -^=5, which it passes through, then goes to the 
other single root at r = (2n + l)7rTo, then returns to the other single 
root (after passing through the double root again), at r = {2n + 2)7rro. 
Neither single root is crossed. The physical time function is; 

t — c~^{bT^/c + Vc + b'^ sm{T^/c)). 

If we start the universe at a big bang with u — 00, we need y — 
0, so 2; = — ^^^^2 • Then we return to a big crunch the next time 
that z returns to its initial value. During that time the quantity T^/c 
increases by 2n, so the full cycle of the universe lasts the finite time 
27r|6|c"5 . The physical value of 6 is — |, giving the lifetime as ttc" 2 = 

^(-2i?)-i. 
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6 The generic quartics: the two periods 

Lastly we need to consider the generic case, where the roots are all distinct. We 
may then write the quartic at first in the following form: 

H{x,y) = {{x-ayy + by^)){{x-cyy + dy^). 

Here a, b, c and d are real and the roots are - = a ± \/^, c ± \/^ and b and d 
must be non-zero. Multiplying out, we find that for the physical quartics, we may 
take: 

c = + 2a), b = ia(l -p),d= ic(l +p) = -^(1 + 2a)(l 

Here a 7^ 0, a 7^ — |, p 7^ 1 andp 7^ —1. Also we have: 

R = -(p(l + 4a) + 2a2 + 8a-l),5= 16a(l + 2a)(l-p + 4a)2. 
9 

lfp= 1 + 4a, then H{x, y) has a factor of x^, contrary to hypothesis. So we have 
p 7^ 1 + 4a. So A is non-zero and has the sign of a(l + 2a). Then we find that the 
four roots are always distinct, except for the following value of p, which must be 
ruled out : 

p= (4a+l)(64a2 + 32a + l). 

This parametrization is symmetrical under (a, p) — > (— a — |, —p), so, without loss 
of generality, we may henceforth assume that a > —\. In particular 2a + 1 > 
and S and A have the sign of a. 

• The four roots are real if and only if a{p — 1) and (1 + 2a) (p + 1) are both 

positive, so if and only if a > and p> 1 (so A > 0) , or — | < a < and 
-1 < p < 1 (so A < 0). 

• The four roots are non-real if and only if the a(p — 1) and (1 -|- 2a) (p -|- 1) 
are both negative, so if and only if a > and p < — 1. In this case A > 0. 

• Exactly two roots are real if and only if a(2a -I- 1) (p^ — 1) < is negative, 
so if and only if a > and — 1 < p < 1 (so A > 0), or — | < a < and 
IpI > 1 (so a < 0). 
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change notation to analyze each of the three cases in turn. 

• If none of the roots is real, then by an appropriate GL(2, M) transformation, 
we can put the quartic in the following standard form: 

H{x,y) = A\x^ + y^){x^ + kY). 

Here < A; < 1 and >1 > 0. Put A;' = Vl - k'^, so < A;' < 1 also. We may 
construct the required transformation, as follows, where we assume that the 
roots are known, for example, using the parametrization given above. Let 
H{x, y) be given initially as follows: 

H{x, y) = {{x - yaf + h''y''){{x - ycf + dV)- 

Here a, b, c and d are real and b, d are positive. Also the four complex roots, 
a±ib and c ± id are assumed to be pairwise distinct. Then we make the 
parametrization: 

[x, y] = [s{a^ + 6^ — bkd — ac) + akd — be, s{a — c) + kd — b]. 

The inverse transformation is: 

x{b — kd) + y{akd — be) 
x{a — c) — y{a'^ + 6^ — bkd — ac) 

Here the real parameter k is required to obey the quadratic equation: 

bd(k^ + 1) - (d^ + c^- 2ca + + b^)k = 0. 

This quadratic equation for the unknown k has discriminant: 

{{{aMb)-{cMd)){{a-ib)-{c-id)){{aMb)-{c-id))){{a-ib)-{cMd))). 

This is always positive, so the two roots for k are real, distinct, positive and 
have product 1, so there is exactly one root for k in the interval (0, 1), which 
gives us the required value of k. 
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Then the parametrization is always well-defined and with respect to this 
parametrization, the equation for r becomes: 



, 9 kds'^ 

Hf^ — 

6d(l + s2)(P + s2)' 



So we have A = Vk~^bd. Then we have 

/k 
— ds — — iEUipticFfis, k~^). 

^^l + s^){k^ + s^) 

For fixed k, the inverse function is: 

s = -iJacobiSN(iA/cT, ■^). 

As r varies over the real line, the ratio x : y runs through all possible values. 
Then the complete integral gives the conformal period T of the universe: the 
amount that r elapses for the universe to return to its original state. We have: 

AkT^k / (i5 = 2A;EllipticK(A:0. 

As an analytic function of r, s is doubly periodic. The second period is 
imaginary: iT', for T' real and is given by times twice the integral the 
integral of along the imaginary axis from one root to another. So here 
we have: 

rk 1 

AkT'^2k / c^g = 4A;EllipticK(A:). 

J-k V(l-s2)(A;2-s2) 

Since as r varies, all ratios y : x are traced out, the integral for the physical 
time function has a pole when a; = and the physical age of the universe is 
infinite. These universes all have A > 0. 
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• If exactly two of the roots are real, then by an appropriate GL(2, R) trans- 
formation, we can put the quartic in the following standard form: 



H{x,y) = A\y' - x'){x' + kY). 

Here k > and A > 0. We may construct the required transformation, as 
follows, where we assume that the roots are known, for example using the 
parametrization of the last section. Let H{x, y) be given initially as follows: 

H{x, y) = {{x - yaf - hY){{x - ycf + dV). 

Here a, b, c and d are real and b, d are positive. Then the four roots, a ± 6 
and c ± id are pairwise distinct. Then we make the parametrization: 

[x, y] = [s{a^ — + bdk — ac) + akd — cb, (a — c)s + kd — b. 

The inverse transformation is: 

x(kd — b) + yibc — akd) 

s ^— 

x{c — a) + y{a? — 6^ + bkd — ac) 

Here the parameter k is required to obey the quadratic equation: 

bd{k'^ - 1) - (6^ - - - + 2ac)k = 0. 

The product of the roots is —1, so both roots are real and exactly one is 
positive, which gives us the required value of k. Then the parametrization is 
always well-defined and with respect to this parametrization, the equation 
for T becomes: 

kds^ 



dr'^ 



M(l-s2)(s2 + py 



So necessarily |s| < 1 and we have A = Vbdk^. The solution can be taken 
to be: 

A;EllipticK(^=i=^) - /cEllipticFfVl^^, -7=^^ 
Akr^EimcF{s,tk~') = ^ ^ 



VT+k 

The inverse solution is: 

s = -iJacohiSN{AkT,ik~^). 
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As T varies, the solution oscillates smoothly between the value s = — 1, 
when we are at the root x : y = {a — b) : 1 and the value s = 1, when we 
are at the root a; : y = (a + b) : 1. The complete integral gives the conformal 
period T of the universe: the amount that r elapses for the universe to return 
to its original state. We have: 

ds Ak 1 

= EUipticKf , ) = 4EllipticK(iA;-^). 

The second period is imaginary iT', for T" real and is given as follows: 

AkT' = 2k r = = , EllipticK( , ^ )■ 

In the case that A > 0, the evolution goes through a; = 0, so the integral for 
the physical time diverges to infinity there and the universe has infinite age. 
In the case that A < 0, the evolution avoids x = 0, so the integral for the 
physical time is finite and the universe has finite age. 

• If all four of the roots are real, then by an appropriate GL(2, M) transforma- 
tion, we can put the quartic in the following standard form: 

H{x,y)^A\y'-x'){x'-kY)- 

Here A; > and ^4 > 0. We may construct the required transformation, as 
follows, where we assume that the roots are known, for example using the 
parametrization of the last section. Let H{x, y) be given initially as follows: 

H{x, y) = {{x - yaf - bY){{x - ycf - d^)- 

Here a, h, c and d are real and h, d are positive. Also the four roots, a ± 6 and 
c± (i are assumed to be pairwise distinct. Finally, without loss of generality, 
we may assume that both the roots a ± 6 are larger than both the roots c ± d, 
so we have: a-\-h>a — h>c-\-d>c — d. 
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Then we make the parametrization: 

[x, y] = [(a^ — — bdk — ac)s — akd — be, (a — c)s — b — dk]. 

The inverse transformation is: 

x{b + kd) — y{akd + be) 
x{a — c) + y{—a^ + 6^ + bkd + ac) 

Here the parameter k is required to obey the quadratic equation: 

f{k) = bd{k^ + 1) + (6^ + - - + 2ac)k = 0. 

The product of the roots is 1. Also we have: 

/(I) = -((a - 6) - (c + dma + b)-{c- d)) < 0. 

So the roots are real and there is exactly one root with A; > 1, which gives 
us the required value for k. Then the parametrization is always well-defined 
and with respect to this parametrization, the equation for r becomes: 

, 9 kds'^ 
dr^ = 



So necessarily 1 < |s| < A; and we have A — Vbdk~^. Then we may take: 



I _f-2 

AkT = EUipticFfW — , Vl - k-^) - EllipticK(Vl - A;-^). 

V 1 — R 

The inverse expression is: 

s-2 = 1 + {k-^ - l)JacobiSN(aA;(T + tq), Vl-k-^f, 

To = (aA;)-^EllipticK(Vl - k'^). 
The real period T is given by the formula: 

AkT = 2 A; / = 2EllipticK(Vl - A;-^). 

V ~ l){k?' — s^) 

The imaginary period is iT' , where we have: 

/■I dt 
AkT' = k / = 2EllipticK(A;-M. 
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There are two possible evolutions, according to whether or not s > 0. 

- If s > 0, as r varies s cycles smoothly between its minimum value 
of 1 and its maximum value, k > 1. At s = 1, we are at the root 
X : y = {a + b) : 1; at s = /c, we are at the root x : y — {c — d) : 1. 
Neither root is crossed during the evolution and the other two roots are 
avoided completely. 

- If s < 0, as r varies s cycles smoothly between its minimum value 
of —k and its maximum value, —1. At s = —1, we are at the root 
X : y = {a — b) : 1; at s — —k, we are at the root x : y = {c + d) : 1. 
Neither root is crossed during the evolution and the other two roots are 
avoided completely. 

In the A < case the universe has finite lifetime. In the A > case, the 
universe has infinite lifetime. 

Finally note a general fact that has emerged in these case by case analyses: that 
near a simple zero of the cosmic quartic, the conformal time parameter is related 
to a spinor variable u which vanishes at that zero, by a formula of the form u — 
c{t — tqY, to lowest order, with c ^ 0. Then the metric is smooth at the zero, 
in terms of the conformal time parameter r and the zero is not crossed. Thus, in 
the generic case that the quartic has no repeated roots, the conformal geometry 
must oscillate perpetually between two adjacent real zeroes of the cosmic quartic, 
if the cosmic quartic has at least one such zero (since then it must have at least two 
such). If the quartic has no real zeros, then the conformal time is unbounded and 
the spinor parameter ha ratios range over the full projective circle. In this case 
we may regard the conformal geometry as being periodic in real conformal time, 
with period given by the integral of y^^TTAdnB over the projective circle. If the 
quartic has a repeated real root, then this root can be crossed. 
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7 Summary and discussion 



We have shown that cosmological metrics can be put in the following standard 
form, using a real projective spinor co-ordinate ha, where ha is equivalent to k-KA 
for any non-zero real k: 

a^TTA 

Here g,^ is a three-geometry of constant curvature k, with = k. Also i?^^ 
is a constant spinor, which is symmetric in the pairs AB and CD and is trace- 
free. Also and are fixed (linearly independent) spinors. The co-ordinate 
T measures conformal time: the vector field dr is the conformal Killing vector 
orthogonal to the surfaces of homogeneity. The physical time variable t is given 
in terms of r by the differential relation: 

dt = u~^dT. 

The equation y^^ y^^ at^cHbd'^ et^ f = 4, with tt^ complex, defines a torus 
in the case that H*^^ is generic. An equivalent but perhaps slightly less elegant 
formulation of the metric is: 

Here g^^ is a flat metric of signature (2, 1) in three-dimensions and H^^ is a constant 
tensor, with H"-'^ = g°-'^Hl symmetric. Also the metric is restricted to the null cone 
in Pa space: g'^^PaPb = and pa is identified with spa, for s > any real number. 

Also U — " , for constant null vectors u"- and v"- (here the connection with the 

PbV 

previous metric is given by using writing pa = hat^b, g"''^ — e^^^e^^^ (where 
e^^ = — e^^), u"- = a"^a^ and = P^P^ . From this viewpoint there are 
obvious questions and possible generalizations: why does U take this form? Why 
is constant? Is the six-dimensional metric of signature (2, 4) given by ds'^, 
but with the constraint that pa be null removed, of some physical interest? Is 
determined from structure at the big bang? At the time of writing we do not know 
the answer to these questions. 
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However it does seem not unlikely that at the big bang, there is some kind of phase 
transition for which the tensor H^^ as an order parameter: in other words that it's 
specific value is not a priori predictable, but it belongs to a class of such parame- 
ters that might be predictable as a whole. Then once the tensor is set, the basic 
disposition of matter, radiation, topology and cosmological constant is imprinted 
for ever on the universe. 

In the generic case we have found that the conformal geometry simply oscillates 
for ever between two roots of the quartic spinor polynomial defined by H^^ (un- 
less there are no real roots, in which case all values of the projective spinor are 
visited repeatedly). This may be relevant for theories of ekpyrotic [7] or cyclic 
universes [8]. The physical geometry then occupies a part of the conformal evo- 
lution. 

Next we have observed that generically, in the complex, the differential form 
defining the time dt — u^^dr has an imaginary period, corresponding to the in- 
tegral along the path in imaginary r going once around the torus defined by H^, 
with the interval of integration in r given by a complete elliptic integral. We have 
shown how to organize the appropriate co-ordinate transformations to put the pe- 
riod integral into a standard form used in elliptic integral theory. Given H^, it is 
possible to compute the period integral in terms of elliptic functions. This imag- 
inary period is unambiguous in the case that A < 0. However in the case that 
A > 0, there is an ambiguity that arises from the contribution of the pole at m = 0. 
One could eliminate the ambiguity if one simply prescribes that the contour is not 
allowed to stray into the the pre-big-bang era: 3?(m) < 0. However it is not clear 
whether this is physically the right thing to do. So at least in the case that A < 0, 
we have a well-defined periodicity in imaginary time, which is directly calculable 
in terms of the ingredients of the cosmic balance equation. Now since the work 
of Hawking on black hole radiation, it has been realized that periodicity in imagi- 
nary time is associated with a mass and via Planck's constant with a temperature. 
Using standard values for the fundamental constants of nature, the formula for 
the mass is 3.2126(10^^) (St) kilograms, where the periodicity Q=(t) in imaginary 
time is measured in seconds. So we have a natural mass scale for such universes, 
in principle measurable from cosmic data. 
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One would like to give an actual value to the parameters it! and S that determine 
our tensor H^. In section three above these parameters are given directly in terms 
of the quantities involved in the cosmic balance equation. However their values 
are very sensitive to the value of the radiation component of the universe, which 
seems to have been neglected in current experimental measurements, because it 
is believed to be small. Instead these measurements have focussed on the matter 
component and on the A term. It should perhaps be emphasized that Einstein's 
equations are inexorable and apply to the entire evolution of the universe. There- 
fore, what is needed, it seems is a global fit for these parameters, taking into 
account the whole history of the universe. Even to determine these parameters 
with an error of less than ten percent would be a terrific achievement. Note that 
there is one combination of these parameters namely SR~^ which is independent 
of the radiation density, so this combination may be easier to pin down. 

Finally it is worth mentioning that conformally flat cosmological models often 
have a nice description in terms of twistor theory [6]. In the present case the 
three-parameter family of world lines of the vector field can be described by 
a three-parameter family of quadrics in projective twistor space. Using standard 
twistor variables and two-component spinors, these quadrics have the equation: 

Here the twistor is Z"' = (a;^,7r^/), the spinor x"^'^' = x^'^' represents the 
particular world-line and the timelike vector ta represents its tangent vector to the 
world-line. After deleting the line at infinity, each quadric (of real dimension four) 
has topology §^ x C. Then the identification in imaginary time corresponds to an 
identification in twistor space of the form (c<j^, tta') = (cc;^ + st^Tr^', tt^/), with s 
real. This wraps up the C-factor giving the quadrics the topology of the product of 
a sphere and a cylinder. If we further identify corresponding to a real periodicity 
in the variable r, then the cylinder factor is wrapped into a torus. 
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